Spin Valve Effect in ZigZag Graphene Nanoribbons by Defect Engineering 
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We report on the possibility for a spin valve effect driven by edge defect engineering of zigzag 
graphene nanoribbons. Based on a mean-field spin unrestricted Hubbard model, electronic band 
structures and conductance profiles are derived, using a self-consistent scheme to include gate- 
induced charge density. The use of an external gate is found to trigger a semiconductor-metal 
transition in clean zigzag graphene nanoribbons, whereas it yields a closure of the spin-split bandgap 
in the presence of Klein edge defects. These features could be exploited to make novel charge and 
spin based switches and field effect devices. 



The recent fabrication of single and few layers 
graphene has sparked considerable expectations in the 
field of all-carbon nanoelectronics [l|, Although 2D 
graphene is a zero gap semiconductor, the possibility to 
pattern graphene nanoribbons (GNRs) with widths of 
few tens of nm, has enabled band gap engineering Q 
and the development of efficient GNRs-based field ef- 
fect transistors However, the further optimization of 
GNRs-based devices or the design of novel device func- 
tionalities demand for more understanding and control of 
not only unavoidable structural disorder and defects, but 
also edge geometry and chemical functionalization. Ex- 
perimental characterization either with Raman studies, 
scanning tunneling microscopy or high-resolution elec- 
tronic transmission microscopy has revealed a large spec- 
trum of topologically different edge profiles, which can 
exhibit either an armchair or a zigzag symmetry , as 
well as a zoology of edge imperfections [6|. In particular, 
individual pending atoms named Klein defects have been 
predicted ^Tj] and more recently observed [8|, (9|] . In prin- 
ciple, such edge disorder is unsuitable for keeping good 
conduction properties of otherwise clean GNRs. As a 
matter of fact, several experimental [lo| and theoretical 
ll| reports evidence the enlargement of the transport 
gap and strong fluctuations of low temperature conduc- 
tance. 

A pioneering work by Son et al. I4I has unveiled half- 



metallic behavior in clean ZGNRs in the presence of a 
transverse electric field, opening perspectives for the de- 
sign of spin-dependent switching devices Several 
recent studies have also reported on the possibility for 
half-metallicity in ZGNRs using chemical doping 
whereas the application of an external magnetic field was 
shown to trigger a transition from parallel to antiparal- 
lel magnetic edges, resulting in giant magnetoresistance 
phenomena and efficient spintronics devices [l5l |. 

In this Letter, we investigate the electronic and con- 
ductance properties of gated zigzag GNRs (ZGNRs) with 
and without the presence of Klein edge defects. The en- 
ergetics of the GNRs is described by a mean-field Hub- 
bard Hamiltonian which has been shown to well repro- 



duce density functional theory results in the spin local 
density approximation [3, [l3|- The gate-injected charge 
density is also self consistently included in the band struc- 
tures and transport calculations. For clean zigzag GNRs, 
the external gate is shown to tune the electronic structure 
from a semiconducting state to a purely metallic state, 
thus switching on current through the device, whereas 
in presence of Klein defects, the spin-split gap is closed, 
resulting in a transition from a pure spin-current to a 
spin-degenerate charge current, close to the Fermi en- 
ergy. 

ZGNRs are very peculiar in the sense that they have a 
magnetic ground state, with each edge having electrons 
aligned ferromagnetically with each other, albeit anti- 
ferromagnetically with respect to the other edge [l^ . 
The ground state is entirely dominated by these mag- 
netic edge states 0, which generate a semiconducting 
band gap [13] ■ The introduction of a Klein edge profile 
(with additional 7r-electron hopping bonds) in the ZGNR, 
is known to produce a flat band over the entire Brillouin 
zone 0, within a tight-binding description. However, 
with the incorporation of electron correlations, a ground 
state with a net magnetization and a spin-split band gap 
emerges 18, 1^. Kusakabe et al. have theoretically dis- 
cussed a Klein edged situation by dihydrogenating one 
of the edges of the ZGNR, validating the occurence of 
a spin-polarized band at the Fermi energy [ll|. Very 
recently, these edges were observed for the first time, al- 
though locally, in graphene using HR-TEM 

Here, it is our interest to study the effect of an external 
gate on the band gaps of these two categories of bipar- 
tite semiconducting systems, one of which has a anti- 
ferromagnetic ground state (ZGNR), and the other fer- 
romagnetic (ZGNR with a Klein edge defect). We model 
GNRs using the spin-unrestricted mean-field Hubbard 
Hamiltonian, 

{ij),a i,<J 

where cj^, Ciu and riic are the creation, annihilation and 
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FIG. 1: (color online) Band structures of the clean ZGNR with (a) Vg — (b) Vg — 1.5 V for both up and down (degenerate) 
spins. Band structures of the ZGNR with one Klein edge when (c) Vg = (d) Vg = 1.5 V for up (red dotted) and down (black 
solid) spins. The schematic shows few unit cells of the clean ZGNR on the left. On the right, three Klein defect edges in an 
otherwise clean ZGNR have been shown for clarity. 



number operators for an electron of spin a in the tt orbital 
of the i*'^ C atom in the ribbon, with t = 2.6 eV and U = 
2.75 eV The unit ceU of the ZGNR has width 20 (40 
atoms) and that of the Klein edged ZGNR has an extra 
C atom at one edge. The source-drain electrodes, for 
simplicity are assumed to consist of the same GNRs. The 
Fermi energy of the non-gated system (Ep) is obtained 
by integrating the charge density up to half- filling. 

One common way to tune the device conductance is to 
use an external metallic gate for monitoring the depletion 
or accumulation of charges in the conducting channel Q . 
The presence of a third (gate) terminal is numerically 
simulated by shifting the on-site energies by eVg where 
Vg is the gate voltage. The modified charge density distri- 
bution up to Ep at every Vg is further self-consistently 
computed. To do so, the excess charge carriers intro- 
duced by the gate are obtained as n{Ep) ~ n{Ep) where 
n{Ep) is the charge density at the Fermi energy in the 
presence of the gate and n{Ep) is that in the absence 
of the gate {niEp) = j'^^p{E)dE, with p{E) the den- 
sity of states). The modified charge density at every Vg 
is then incoporated into the Hamiltonian and solved self- 
consistently, to obtain the new band-stucture of the gated 
GNR, with the associated Fermi level {Ep). 

The coherent transport through the system is then 



calculated, using the well-known Landauer's formalism. 
The retarded Green function of the device is computed 
as G = [E — H — Tjl — E_r)^^ where and Y^b. are 
the self-energies of the GNR electrodes, calculated as 
S = rgsT^ where r contains the device-source/drain in- 
teractions and is the surface Green's function of the 
electrodes, calculated using standard recursive Green's 
function techniques 11]. The conductance of the de- 
vice is finally obtained as Gd{E) = ^T{E) where T(E), 
the transmission probability at energy E is given by 
Tr[^LG^i^Gt], Tr is the Trace and Tl r ^ ii^L r ~ 

Electronic band structures of ZGNR with or without 
a Klein edge, and as a function of applied external gates 
can be seen in Fig 1. First, as observed in Fig la, when 
Vg = 0, the ground state is characterized by a band 
gap, resulting from its oppositely spin-polarized edges. 
With the application of a positive gate voltage channel- 
ing holes into the system, the electron density at the 
edges (2/3 < \k\ < 1) gradually diminishes, which results 
in a completely non-magnetic, metallic system. A closing 
of the band gap is seen in Fig lb for Vg = 1.5 V. The 
effect of gate voltage on the low-energy conductance can 
be appreciated in Fig 2. Given that the GNRs are disor- 
der free, the conductance remain quantized and is fixed 
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FIG. 2: Spin-degenerate (up and down) conductance of the 
ZGNR with respect to the incident energy of the electrons (E) 
scaled by Ef, for (a) K, = (b) l/g = 1.5 V. 
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FIG. 3: (color online) (a) Variation of the bandgap over the 
entire Brillouin zone with Vg (b) The total number of up (red 
dashed) and down electrons (black solid) with respect to Vg. 
Dotted lines are marked ai Vg = and Vg = 1.5 V for clarity. 



by the number of available conduction channels at the 
relevant Fermi level. At injection energies close to Ep, 
a current switch is driven by the gate voltage increase, 
indicating a transistor type behavior. 

In contrast, the situation for the ZGNR with one Klein 
edge shows some marked differences. First by looking at 
the band structure at zero gate voltage (Fig.lc), it is 
clear that the highest occupied band is made up entirely 
of one spin (dotted) and the lowest unoccupied band, of 
the other (solid). The ground state of the system shows 
some ferromagnetism, with both edges displaying high 
electron density of the same (majority) spin. This is ex- 
plained by the finite sublattice imbalance which favors 
the appearance of midgap states and magnetic proper- 
tics [16, 17, 20]. 

When a positive gate voltage is applied, the major- 
ity spin (up) band shifts upward, reducing the electron 
(up spin) density, as seen in Fig 3b, which eventually 
meets the minority spin (down) density to make the sys- 
tem completely non-magnetic and metallic. This is also 
characterized in Fig 3a by the slow reduction of the spin- 
split gap, which eventually disappears beyond Vg = 1.5 
V. The bandstucture of the GNR at Vg = 1.5 V, shown 
in Fig Id, clearly exhibits spin-degenerate bands. When 
a negative gate bias injecting excess electrons into the 
device is applied, a similar situation follows. The down 
spin band gets pushed into the occupied states, dimin- 
ishing the total magnetization and hence the spin-split 
gap. Here, the down spin density increases to meet its 
up spin counterpart as seen in Fig 3b. 

In Fig. 4 at Vg = 0, the transmission probability of the 
majority spin close to Ep is very high, as it is the only 
band available in the energy range. As the gate voltage 
is increased, the system transmits only the majority spin 
until ss 1.0 V, after which it turns non-magnetic with 
both spins having identical transmissions. This indicates 



that it is possible to tune the current in the system from 
one which is completely spin-polarized to an unpolarized 
charge current just with the help of an external gate. 

The origin for such transitions from semi-conducting to 
metallic in the gated ZGNRs and from a spin-polarized 
to a spin-degenerate response in the Klein edged ZGNRs 
can be rationalized as follows. The ground state of both 
systems are entirely composed of magnetic edge states, 
which produce their respective band gaps. An exter- 
nal gate voltage does break the spin polarization of the 
system due to the change of the charge number at the 
edges. As a consequence, a controlled switch from a 
spin-splitting gap to metallicity is observed. One notes 
that the application of an external magnetic field in clean 
ZGNR a lso yields a similar effect and giant magentore- 
sistance 



17|. 



In summary, we have studied the effect of gate volt- 
age on the band-structure and conductance profiles of 
a ZGNR and a Klein edged ZGNR, based on the mean- 
field spin-unrestricted Hubbard model. The gate-injected 
charge density was obtained sclf-consistently, and con- 
ductance calculations were performed using the Lan- 
dauer's formalism. We have found that an external gate 
could tune the spin-induced band-gap in clean ZGNRs or 
in Klein edged ZGNR. In the former case, this leads to a 
switching-on of the charge current, whereas in the latter 
case, a transition from a pure spin current to a com- 
pletely unpolarized charge current is achieved. These 
features could eventually help in designing GNR based 
charge and spin switches and FETs in future all-carbon 
circuits. 
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FIG. 4: (color online) Conductance of the Klein edged ZGNR 
with respect to the incident energy of the electrons (E) scaled 
by Ef, for (a) Vg = (b) Vg = 1.5 V. Dashed (red) lines refer 
to the spin-up conductance and solid (black) lines refer to the 
spin-down conductance. 
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